Rheology of dilute suspensions of vesicles and red blood cells 
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We present rheology experiments on dilute solutions of vesicles and red blood cells (RBC). Varying 
the viscosity ratio A between internal and external fluids, the microscopic dynamics of suspended 
objects can be qualitatively changed from tank-treading (tt) to tumbling (tb). We find that in the 
tt regime the viscosity rj, decreases when A increases, in contrast with droplet emulsions and elastic 
capsule theories which are sometimes invoked to model RBC dynamics. At a critical A (close to the 
tt-tb transition) rj exhibits a minimum before it increases in the tb regime. This is consistent with a 
recent theory for vesicles. This points to the nontrivial fact that the cytoskeleton in RBC does not 
alter the qualitative evolution of and that, as far as rheology is concerned, vesicle models might 
be a better description. 

PACS numbers: 87.16.Dg 47.57.-s 87.19. Tt 47.63.-b 
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Complex fluids are ubiquitous in nature, and their 
study embraces a large spectrum of disciplines: physics, 
biology, chemistry, engineering... Prominent examples in 
biology are blood and cartilage. Complex fluids are also 
essential in several domestic applications: cosmetic liq- 
uids (e.g. shampoo), emulsions (e.g. mayonnaise), sus- 
pensions (e.g. clay ), etc... 

Ordinary fluids and elastic solids obey universal laws, 
namely Navier-Stokes and Hooke-Lame laws. Contrari- 
wise, complex fluids still continue to escape a univer- 
sal description (if any). The difflculty lies in the in- 
timate coupling between the microscopic dynamics and 
the macroscopic flow. For ordinary fluids, like water, the 
molecular time scale (molecular rotation, and vibration) 
is in the range of 10~^^ s, which is much smaller than 
macroscopic flow scales, so that an adiabatic elimination 
of molecular modes in favor of hydrodynamical ones is 
legitimate. In contrast, for a complex fluid like blood, 
red blood cells (RBC) move and deform over time scales 
comparable to the global flow, so that the macroscopic 
constitutive law should carry information on the dynam- 
ics of the suspended entities, even if the law is averaged. 

The understanding of some complex features of the 
rheology of vesicle and RBC suspensions has recently 
progressed thanks to a coupling between experiments and 
modeling. Our experimental study focuses on the rheol- 
ogy of dilute suspensions of vesicles and RBC. The study 
of vesicles under nonequilibrium conditions keeps receiv- 
ing an increasing interest. They constitute a relatively 
simple system (albeit complex enough at the absolute 
level) that is believed to capture some features of RBC. 

We choose to study a dilute suspension for the follow- 
ing reasons, (i) This is a simple starting point which 
avoids additional complexity like hydrodynamic interac- 
tions between objects or rouleaux formation for RBC's 
(ii) As for other complex fluids (e.g. dilute polymer solu- 



tions where the celebrated Oldroyd B model constitutes 
a basic reference[l]) this offers a favorable terrain for con- 
frontation between experiments and modeling based on 
a full microscopic description of the suspended entities. 
By focusing simultaneously on vesicles and RBC, we can 
draw information regarding rheology, about similarities 
and differences of the two systems. 

A key parameter in the present study is the viscosity 
contrast A = rjin/rjout, where rjin and rjout are the vis- 
cosities of the internal (hemoglobin for RBC, sugar -f 
dextran solution for vesicles) and the external solution 
(e.g. buffer -I- dextran, water -I- dextran). 

Vesicles are known to exhibit three kinds of motions: 
(i) tank-treading in that the vesicle orients steadily 
itself with a certain angle with respect to the flow direc- 
tion (while its fluid membrane makes a tank-treadingmo- 
tion), (ii) tumbling when A exceeds a critical valueJ^H], 
and (iii) vacillating-breathing (vb) mode where the vesicle 
long axis oscillates (or vacillates) around the flow direc- 
tion, whereas its shapes undergoes a breathing motion( 
see [1, [1,0, 11)0 for theory, and [lol, [ll| for experiments.). 
An interesting question immediately arises: do the micro- 
scopic dynamics {tt, tb or vb) and the transitions between 
them have a signature on the macroscopic level of rheol- 
ogy? 

A major result reported here is the fact that the experi- 
mental rheology of RBC suspensions shows the same gen- 
eral trend as the theory for vesicles: the effective viscosity 
exhibits a minimum in the vicinity of the tt-tb transition. 
On the one hand, this reveals a qualitative change due to 
the link between microscopic and macroscopic dynamics. 
On the other hand, this suggests that the cytoskeleton 
does not induce a significant qualitative change on rhe- 
ology. This behavior is consistent with recent theoretical 
studies on vesicle suspensions [§, [12], and contrasts with 
other theories sometimes invoked for RBC dynamics, like 
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droplet [13| and capsule [lj| theories, where the viscosity 
is predicted to increase with A. 

Theory: Einstein [l5| provided the famous expression 
of the effective viscosity of a dilute suspension of rigid 
spherical particles. The intrinsic viscosity is given by 



(1) 



where 779 is the solvent viscosity and the volume fraction 
of particles. Later, Tavlor[13j provided the analogous 
expression for an emulsion: 



W\ = 



5A/2 + 1 
A + 1 



(2) 



A step further consists in adding a membrane around 
the drop. This membrane can either be an elastic solid 
(capsule) for which expressions for the viscosity of the 
suspension have been suggested or a fluid lipid bi- 
layer (vesicle). Recently an expression has been derived 
in the tt regime for quasi-spherical vesicles [B^ 
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(the subscript stands for tank-treading motion), where A 
is the excess area relative to a sphere A — [A — A-kt"^) / 
where r is the sphere radius. Several noticeable differ- 
ences with droplets and capsules can be mentioned: (i) 
the viscosity decreases with A for vesicles, while the con- 
trary is found for droplets and capsules, (ii) When A is 
large expression ([3]) does not tend to 5/2, as does ©. 
Expression [3] is valid in the tank-treading regime only 
0. At low shear rateQ there is a direct bifurcation from 
tt to tb when increasing A. In this limit, one can as- 
sume a shape-preserving motion. Following the general 
expression for the instantaneous viscosity [12], [77] can be 
derived in the tb regime (technical details will be given 
elsewhere) : 
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with h = 60^27r/15/(23A + 32). The tumbhng domain 
corresponds to Ah^ < A (the opposite limit is the domain 
of tt)[5|. Figure [3] shows the behavior of [rj] for vesicles in 
both the tt and tb regimes. At the bifurcation, one has 
a cusp singularity, with a linear behavior on the tt side, 
and a square root singularity on the tb side. 

Materials and methods: Viscosities of diluted vesicu- 
lar and red blood cell suspensions were measured as a 
function of the viscosity ratio A between the inner and 
outer fluids (for the vesicle or erythrocyte membrane). 
Chemicals were purchased from Sigma-Aldrich (France) 
and used without any further purification. Microscopic 
observations were performed on an Olympus IX 71 in- 
verted microscope in phase contrast (GUVs and RBC) 



or bright field (RBC) modes. An Eppendorf 5804 cen- 
trifuge was used during sample preparation. 

Erythrocyte samples were provided by the CHU (Cen- 
tre Hospitaller Universitaire) of Grenoble (France) from 
haematologically healthy donors. Phosphate buffer saline 
(PBS, containing 0.01 M phosphate buffer (10.1 mM 
Na2HP04 and 1.8 mM KH2PO4), 0.0027 M KCl and 
0.137 M NaCl) with pH 7.4 and osmolarity 290 ± 10 
mOsm/kg or PBS dextran solutions have been used as 
suspending media. Blood samples were gently centrifu- 
gated at 1200 g for 5 min and the buffy coat was re- 
moved. Erythrocytes were washed several times that way 
with PBS and the hematocrit was measured via a stan- 
dard automated analyser. The healthy biconcave shape 
of erythrocytes in the studied suspensions has been sys- 
tematically observed in phase contrast or bright field mi- 
croscopy. The cytoplasmic viscosity of RBC is deter- 
mined on the basis of the mean corpuscular haemoglobin 
concentration (MCHC) [2l|. At 22°C it was estimated 
to be of the order of 20 mPa.s [l^]- The final samples 
were prepared by diluting starting suspensions into a so- 
lution of dextran (D4751 and D1037 from Sigma, with 
respective molecular weights of and 64 000-76 000 and 
425 000-575 000) in PBS buffer or pure buffer. Because 
of the aggregation effect of dextran on RBC 18, 3, 2o| . 
only dextran polymers with concentrations higher than 
7% w:v (7 g / 100 ml) were dissolved in the buffer in 
order to modify A. Thus, adding up to 25% of dextran 
in the suspending medium, A was varied in the range of 
0.12-3.41 and 18.37-20 without dextran. 

Vesicles were prepared using the electroformation 
method [l3| with electroformation chambers designed to 
maximize the vesicle yield and concentration. We used 
DOPC, or (dioleoyl phosphatidylcholine, P-6354 from 
Sigma), and the electroformation solution is a sucrose 
solution with a concentration between 200 and 700 mM 
in a 1:4 glycerol- water (w:w) mixture. Samples were di- 
luted (about 7 times) in an outer medium made of glucose 
and dextran in a 1:4 glycerol- water mixture. This outer 
medium is slightly hyper osmotic in order to deflate vesi- 
cles, has a different viscosity in order to vary A and a 
slightly smaller density. Diluted samples were then cen- 
trifugated (at 500-3000 rpm, corresponding to accelera- 
tions of 40-1500 g, which gives rise to sedimentation ve- 
locities and shear rates that are safe for vesicles) . A range 
of A between 0.2 and 1.2 was explored. In principle, it is 
possible to produce vesicles with a more viscous interior 
solution, and therefore reach higher values of A by adding 
dextran in the electroformation solution. However elec- 
troformation in dextran solution is much less efficient in 
two ways: vesicles are smaller and less concentrated and 
the dextran content is extremely heterogeneous, adding 
polydispersity in A in addition to the existing polydis- 
persity in size and excess area of vesicle samples. The 
corresponding data points have therefore been discarded 
and are not presented in this paper. As a general rule. 
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the preparation of large and concentrated vesicle samples 
is lengthy and difficult. 

Viscosity measurements were made using two rheome- 
ters: a stress-controlled Bohlin Gemini 150 rheometer 
(Malvern Instruments, Germany) with a cone-plate ge- 
ometry (60 mm diameter, 2° angle) and a LS30 low- 
shear rheometer (Contraves, Germany) with a cylinder- 
Couette geometry used for cases where the solvent viscos- 
ity is below 7 mPa.s. With both rheometers, the viscosity 
was measured across a range of shear rates between 1 and 
a few hundred (for the LS rheometer: up to the up- 
per limit of torque measurement). All viscosity measure- 
ments with GUV and RBC suspensions were performed 
at the constant temperature of 22°C. 

For RBC samples, solvent viscosity was measured after 
gently centrifugating until a completely clear supernatant 
is obtained. The absence of any suspended RBC has 
been systematically controlled microscopically. For vesi- 
cle samples, an equivalent outer medium was prepared 
following the same dilution procedure as the sample with 
pure interior solution. 

Volume fractions of vesicles and RBC were estimated 
using different techniques. For RBC's, the final hemat- 
ocrit was checked by counting cells in a sample placed in 
a 100 /im deep PDMS observation chamber. During the 
injection of RBC samples in that chamber, we could also 
check that the two regimes of red blood cell dynamics 
{tt and tb) were microscopically observed for A < 2 and 

A > 2 [H,!!!. 

For vesicles, after the viscosity measurement, the sam- 
ple was placed in the observation chamber (after con- 
trolled dilution to avoid vesicle overlapping on the im- 
age), and a statistical measurement of the vesicle pop- 
ulation (vesicle diameter, number) was done using NIH 
Image or ImageJ software. Although this technique re- 
quires a lot of manipulation and image processing, we 
found that this was the most accurate way of determin- 
ing vesicle volume fraction, in contrast with centrifuga- 
tion for instance. 

Results: Viscosity measurements were made in a large 
range of shear rates in order to detect possible viscoelas- 
tic or aggregation effects and determine the range of shear 
rates were measurements are accurate and correspond to 
the dilute suspension limit. For every sample, we re- 
tained the average value of viscosity measured in in the 
plateau at higher shear rates (between 10 and 100 s~^ 
in the example of figure [1]) . The fluctuations or differ- 
ent behavior at low shear rate are either due to a lack 
of precision of the rheometer in the corresponding torque 
range or to aggregation (rouleaux formation) for RBC's. 

For RBC's in PBS buffer, several measurements were 
made at different values of the concentration (hematocrit 
or volume fraction) in order to check the linearity of the 
relationship between viscosity and concentration, and the 
applicability of dilute suspension theories. The results 
are shown in Fig. [2]where up to experimental precision no 
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FIG. 1: Viscosity rj vs shear rate -y for suspensions of red blood 
cells in PBS buffer (A ~ 20) for various volume fractions ( (♦) 
3%, (A) 5%, (.) 7%, (□) 10%). 




FIG. 2: Viscosity (relative to solvent viscosity) vs volume frac- 
tion for suspensions of red blood cells in PBS buffer (A ~ 20), 
showing a linear dependency across nearly the whole range of 



nonlinear behaviour can be detected for concentrations 
up to 10 %. 

Figure [3] shows the measured intrinsic viscosities [rj] as 
a function of the viscosity ratio A for RBC's with </) ~ 5%, 
and vesicles with (j) ranging from 1.87 to 10.5 %. Remark- 
ably, the rheology of RBC suspensions follows the qual- 
itative trends predicted by vesicle theory: a measurable 
decrease of [77] in the tt regime, a minimum near the bifur- 
cation to tb, and a sharp increase in the tb regime when A 
is increased. In Fig. [Sj the viscosity predicted by eqs. [3] 
and [4] with A = 1.5 is shown to reproduce rather well the 
measured viscosity for RBC's, although their excess area 
is closer to A = 4. One should keep in mind however that 
eqs. [3] and [4] approximate solutions for a nearly spher- 
ical vesicle model, that does not take into account the 
capillary number (ratio of viscous forces and membrane 
bending forces). Droplet and capsule models fail to pre- 
dict the correct qualitative behavior in the tt regime and 
do not even approach the measured values of [77] for high 
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FIG. 3: Intrinsic viscosity [77] vs viscosity ratio A for RBC's 
with ~ 5% (•), vesicles (o), and various models (dotted 
line: drops, dashed line: capsules with e = 0.1 [l^, plain line: 
vesicles with A = 0.5 and A = 1.5). 



viscosity ratio, where they saturate at [rf\ — 2.5 like sus- 
pensions of sohd spheres. For vesicles in the tt regime, 
no clear tendency can be observed due to the disper- 
sion and lack of accuracy of data points. Several factors 
are responsible for this: (i) there is a polydispersity of 
vesicle's excess area in the samples, (ii) vesicle's sizes are 
polydisperse, leading to different capillary numbers in the 
same sample, (iii) for most samples, the volume fraction 
is rather small (a few percent), leading to small absolute 
viscosity variations and amplifying errors when comput- 
ing the effective viscosity. However, Fig. [3] reveals that 
the measured values of [77] are compatible with vesicle 
theory [5] for A = 0.5, a value consistent with measured 
excess areas of vesicles, and above the values obtained 
for RBC's, which have a larger excess area. 

Remarkably, the rheology of RBC suspension can be 
described by vesicle theory even though it is a simpler ob- 
ject, which seems to capture the essence of RBC dynam- 
ics under flow, as far as rheology is concerned. A clear 
transition between the tt and th regimes is observed, with 
a minimum effective viscosity around the bifurcation. 
The data obtained for vesicle suspensions is compatible 
with the model and confirms the predicted dependency 
on excess area. For a more detailed and quantitative ex- 
perimental study of the rheology of vesicle suspensions, 
an important step forward in the production of controlled 
samples is a prerequisite. A decisive improvement would 
be the possibility to drastically reduce the polydispersity 
of vesicle size, excess area and viscous content. 
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